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Molecular phylogeny of Magnoliaceae based on plastid DNA sequences with special emphasis on 

some species from continental Southeast Asia 


HIROSHI AZUMA*, PIYA CHALERMGLIN** & HANS PETER NOOTEBOOM*** 


ABSTRACT. A molecular phylogenetic tree of Magnoliaceae was constructed based on four cpDNA regions ( trnK intron including 
matK, trnH-psbA, atpB-rbcL , and ndhF; total ca 5800 bp) to reveal the phylogenetic positions of some Magnolia species from 
Peninsular Malaysia, Thailand, and Vietnam. A total of 100 accessions, including 75 Magnolia, two Liriodendron species, three 
putative hybrids (including M. x alba), and two outgroups ( Degeneria and Galbulimima ), were included, of which 20 species plus 
three putative hybrids (total 34 accessions) were collected in continental Southeast Asia. Phylogenetic analyses (maximum parsimony 
and maximum likelihood) indicated that most of these newly examined species were placed in expected taxonomic groups (e.g., M. 
villosa in subsection Blumiana; M. garrettii, M. hookeri, and M. utilis in section Manglietia; M. thailandica in section Kmeria', M. 
mediocris, M. scortechinii, and M. sirindhorniae in subsection Michelia; M. macklottii var. beccariana in subsection Maingola; and 
M. pahangensis in subsection Aromadendron). Circumscription of the taxonomic groups matched the phylogenetic clades except in 
subsections Gwillimia and Blumiana. Phylogenetic relationships among major large clades were still unresolved and interspecific 
relationships within some large clades (e.g., subsection Michelia and section Manglietia ) were also not suggested. 
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INTRODUCTION 

Magnoliaceae comprise ca. 220 species of 
trees or shrubs distributed in temperate and tropical 
regions from the Himalayas to East Asia and 
Southeast Asia and the Americas, and is divided 
into two subfamilies, Liriodendroideae and 
Magnolioideae (Nooteboom, 1993; 2000; Figlar & 
Nooteboom, 2004). The subfamily Magnolioideae 
had formerly been divided into six (or more) genera 
(Magnolia L., Manglietia Blume, Michelia L., 
Elmerrillia Dandy, Kmeria Dandy, Pachylarnax 
Dandy: Chen & Nooteboom, 1993; Nooteboom, 
1993; Frodin & Govaerts, 1996), but several mor¬ 
phological and molecular studies have shown all 
genera to be embedded in Magnolia, thus the family 
now consists of just two genera, Liriodendron L. 
and Magnolia (Figlar & Nooteboom, 2004). 

Molecular phylogenetic analysis based on 
DNA sequence data has been found to be very useful 
for understanding phylogenetic relationships 


among species as well as evolutionary trends of 
morphological and other taxonomic characters 
(e.g., chromosome number and chemical compounds) 
including geographical distribution. Thus, it has 
become a powerful tool for construction and/or 
evaluation of taxonomic systems (Judd et ah, 
2008). Many molecular phylogenetic studies were 
conducted on Magnoliaceae by using plastid (chlo¬ 
roplast) (Azuma et ah, 1999; 2001; 2004; Shi et ah, 
2000; Kim et ah, 2001; Suh & Kim, 2001; Wang et 
ah, 2006) or nuclear DNA sequences (Nie et ah, 
2008). Several early studies gave many insights 
into the recent classification of Magnoliaceae 
(Figlar & Nooteboom, 2004; Nooteboom & 
Chalermglin, 2009). 

Although these molecular studies commonly 
recognized some well-supported clades, especially 
at terminal nodes of the trees, we have not yet 
attained a fully resolved phylogenetic tree of 
Magnoliaceae because of a limited number of 
species analysed and low resolution among major 
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clades in the tree. For example, Kim et al. (2001) 
constructed a phylogenetic tree of 99 species by 
using only ndhF sequences (ca 2200 bp). This 
resulted in a tree with little or no resolution among 
the major clades. By contrast, Suh & Kim (2001) 
analyzed a combined plastid DNA sequence data 
set of ndhF , rbcL, trnK intron including matK, 
trnH-psbA spacer region, atpB-rbcL spacer regions, 
trnL intron, trnL—trnF spacer region and ORF350 
regions (total ca. 8700bp). The resulting tree did 
slightly improve resolution among the major 
groups, but the tree consisted of only 48 taxa 
(mostly temperate species) of Magnoliaceae. 
Therefore, it is important to continue additional 
molecular work to fully understand the phylogenetic 
systematics of Magnoliaceae, especially in regard 
to tropical species most of which have not yet been 
included in the molecular phylogenetic analyses. 
Moreover, new species are still being discovered in 
the tropics, for example M. sirindhorniae Noot. & 
Chalermglin andM thailandica Noot. & Chalermglin, 
both of which have been discovered in Thailand 
during the past decade (Nooteboom & Chalermglin, 
2000; 2002; see also Figlar, 2011). Previous molecular 
studies dealt with only a few species widely dis¬ 
tributed or cultivated in Southeast Asia [e.g., M. 
champaca (L.) Baill. ex Pierre, M. insignis Wall., 
and M. liliifera (L.) Baill.], thus phylogenetic 
relationships among most Magnolia species found 
in this area are poorly understood. 

The purpose of this study is to construct a 
molecular phylogenetic tree of Magnoliaceae and 
to know the phylogenetic positions of especially 
tropical Asian species. In total, 75 species plus 
three putative hybrids (a total of 98 accessions) 
including 20 newly analyzed species collected 
from Peninsular Malaysia, Thailand, and Vietnam 
were analyzed using plastid DNA sequence data of 
four regions ( trnK vul J intron including matK gene, 
ndhF gene, tmH GVG —psbA and atpB-rbcL spacer 
regions; total ca 5800 bp). We also tried to evaluate 
the taxonomic system of Magnoliaceae by Figlar 
and Nooteboom (2004) in which the genus Magnolia 
is divided into three subgenera and twelve sections 
(thirteen subsections in five sections; Table 1). In 
this study we included representatives of all 
subgenera, sections, and subsections except for 
subsection Dugandiodendron (Lozano) Figlar & 
Noot. (tropical American group) in section Talauma 
Baill. 


MATERIALS AND METHODS 

A list of species (samples) involved in this 
study is shown in Table 1 giving information or 
references, and DDBJ/EMBL/GenBank accession 
numbers of sequences. Most sequences were 
obtained from DDBJ (DNA Data Bank of Japan), 
but sequences of 20 species plus three putative 
hybrids collected in Peninsular Malaysia, Thailand, 
and Vietnam (total 34 samples) were newly deter¬ 
mined in this study. We also determined ndhF 
sequences of 42 samples (taxa) which have been 
used in previous molecular studies(Azuma et al., 
1999; 2001) in which the other three regions ( trnK 
intron, trnH-psbA, and atpB-rbcL ) had already 
been sequenced (Table 1). We also retrieved or newly 
sequenced corresponding regions of Degeneria 
vitiensis I.W.Bailey & A.C.Sm. or D. roseiflora 
John M.Mill. (Degeneriaceae) and Galbulimima 
belgraveana (F.Muell.) Sprague (Himantandraceae), 
which were added to the data matrix as outgroup 
taxa, because they are members in Magnoliales and 
seem to be more closely related to Magnoliaceae 
than Annonaceae and Myristicaceae (Soltis et al., 
2000; Hilu et al., 2003; Qiu et al., 2010). Finally, 
DNA sequence data of a total of 100 accessions 
from 77 species and three putative hybrids including 
M. x alba were compiled and analysed (Table 1). 

Sequences of the following four regions of 
plastid DNA were retrieved or sequenced; the 
tr/7K uuu intron including the matK gene (ca 2500 
bp), the ndhF gene (ca 2090 bp), and trnFF VG - 
psbA (ca 430 bp) adn atpB-rbcL (ca 800 bp) spacer 
regions. These regions have often been used in 
phylogenetic analyses of Magnoliaceae (e.g., 
Azuma et al., 1999; 2001; 2004; Kim et al., 2001). 

Total DNA was extracted from silica-gel 
dried leaves or herbarium specimens using the 
modified CTAB method by Doyle and Doyle 
(1987). The PCR mixture (20 pL) contained 1 pL 
of template DNA, 2 pL of dNTPs (2.5 mM each), 1 
pL of each primer (10 pM), 2 pL of lOx Taq buffer 
(containing 20 mM MgC12), 0.5 U of Taq polymerase 
(TaKaRa Ex Taq, Takara Bio Inc., Japan). The PCR 
was performed with a GeneAmp PCR System 2700 
(Applied Biosystems Japan Ltd., Japan) starting at 
94°C (5 min), followed by 35 cycles of denaturation 
at 94°C (30 sec), annealing at 50°C (30 sec), and 
extension at 72°C (30 sec or 1 min), and a final 
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extension at 72°C (7 min). After checking for a 
single band by electrophoresis on an 1% agarose 
TAE gel stained with ethidium bromide, the PCR 
products were purified by enzymatic treatment. 
That is, 2 pL of mixture containing 0.3 Us of 
Exonuclease I (Takara Bio Inc., Japan) and Calf 
Intestine Alkaline Phosphatase (TOYOBO Co., Ltd., 
Japan) was added to each PCR tube to degrade 
remaining primers and to dephosphorylate remaining 
dNTPs. The tubes were heated at 37°C for 30 min 
then at 80°C for 15 min. Direct sequencing of both 
DNA strands was conducted on an ABI 3100 
Genetic Analyzer (Applied Biosystems Japan Ltd., 
Japan) using a BigDye Terminator version 3.1 Cyclic 
Sequencing Ready Reaction Kit (Applied Biosystems 
Japan Ltd., Japan) following the manufacturer’s 
protocol. Primers used for amplification and/or 
sequencing are listed in Table 2 and the relative 
positions are illustrated in Fig. 1. 

Alignment of sequence data was manually 
carried out in MEGA 4 (Tamura et al., 2007). Four 
regions were combined as a single data matrix. 
Three and one short inversions were found in the 
matK region (all 33 species in subgenera Yulania 
and Gyrtopodium, and section Kmeria showed - 
TCTATT- instead of-AATAGA-; 20 species, i.e., 
M. dawsoniana, all 10 species in section Magnolia, 
all 8 species in section Gwillimia, and L. chinense 
showed -GAA- instead of -TTC-; and 2 species, 
M. nitida and M. sinica, showed -GA- instead of 
-TC-) and in the tmH-psbA spacer region [7 species, 
i.e., M. figo, M. vrieseana (syn. M. ovalis), M. 
acuminata, M. nitida, M. fraseri ssp. fraseri and 
ssp. pyramidata, M. delavayi, and L. tulipifera, 
showed -TTCTAT- instead of -ATAGAA-], 
respectively (Azuma et al., 1999). These inversions 
were excluded from the data matrix. Samples 
which showed completely identical sequence were 
represented by one OTU (operational taxonomic 
unit) in phylogenetic analysis. Indels were treated 
as missing data. 

Maximum parsimony analysis (MP) was 
conducted in PAUP* 4.0bl0 (Swofford, 2003). For 
the heuristic search we used a random stepwise 
addition sequence with TBR (tree bisection and 
reconnection) branch swapping and 10 replication 
of the random addition sequence. One thousand 
bootstrap replications were performed using the 
same setting as in the heuristic search. Maximum 


likelihood analysis (ML) was conducted in MEGA 
5.05 (Tamura et al., 2011). The best substitution 
model was selected based on BIC (Bayesian 
Information Criterion) values of 24 possible 
combinations of substitution models and rate 
description analyzed in MEGA. The GTR+G 
model, the lowest BIC value, was choosen for our 
data set. For ML analysis, gaps were partialy deleted 
by 95% “site coverage cutoff’ option. NNI 
(Nearest-Neighbor-Interchange) was used as ML 
heuristic method and the initial tree was made 
automatically. One hundred bootstrap replications 
were performed using the same setting. 

RESULTS AND DISCUSSION 

A combined data matrix of four plastid 
regions was composed of 5960 characters after 
alignment (total 85 ingroup and 2 outgroup OTUs). 
There are 785 (13.2%) variable sites among all 
OTUs, of which 374 (6.3%) were parsimony- 
informative, and 420 (7.0%) variable sites among 
Magnolia species, of which 223 (3.7%) were 
parsimony-informative. The MP analysis generated 
334 equally most parsimonious trees with a tree 
length 999; the consistency index (Cl) and the Cl 
excluding uninformative characters were 0.8649 
and 0.7648, respectively, and the retention index 
(RI) was 0.8985. The strict consensus tree of 334 
MP trees is shown in Fig. 2 with subgenera, 
sections and subsections of the Magnoliaceae 
according to Figlar and Nooteboom (2004). 
Maximum likelihood tree (InL = -14976.16) is 
shown in Fig. 3. There is no topological contradiction 
between the MP and ML trees. 

The Magnoliaceae were first divided into two 
clades which correspond to the two subfamilies, 
i.e., Liriodendroideae (only Liriodendron) and 
Magnolioideae (only Magnolia) (Fig. 2). In the 
clade of the Magnolioideae, three well supported 
clades [clade B-D, >98% bootstrap percentage] 
and one moderately supported clade (clade A, 79% 
in MP tree and 80% in ML tree) were recognized. 
The relationships among these clades were not 
suggested in the MP tree(Fig. 2), but two tropical 
American clades (C and D) formed a weakly supported 
clade (60%) in the ML tree (Fig. 3). Clade A is a 
relatively large clade and contains taxa from sub¬ 
genera Magnolia [except for sects. Gwillimia DC. 


MOLECULAR PHYLOGENY OF MAGNOLIACEAE BASED ON PLASTID DNA SEQUENCES WITH SPECIAL EMPHASIS ON 
SOME SPECIES FROM CONTINENTAL SOUTHEAST ASIA (H. AZUMA, P. CHALERMGLIN & H. P. NOOTEBOOM) 


151 


(tropical Asia; clade B) and Talauma (tropical 
America; clades C and D)], Yulania (Spach) Rchb., 
and Gynopodium Figlar & Noot., indicating that 
subgenus Magnolia is at least paraphyletic. All 
temperate and some subtropical to tropical species 
were found in this clade. 

Monophyletic relationship of members in 
sect. Gwillimia (clade B), sect. Talauma subsect. 
Splendentes (Dandy ex A.Vasquez) Figlar & Noot. 
(clade C), and sect. Talauma subsect. Talauma 
(clade D) were suggested, respectively. Section 
Gwillimia is taxonomically divided into two sub¬ 
sections, Gwillimia and Blumiana (Blume) Figlar 
& Noot. (Figlar & Nooteboom, 2004). Flowever, 
neither representatives of subsect. Gwillimia nor 
representatives of subsect. Blumiana form well 
supported clades. That is, M. pterocarpa Roxb. and 
M. henryi Dunn belong to subsect. Gwillimia with 
M. coco (Lour.) DC. andM delavayi Franchet., but 
the former two species seem to be more closely 
related to M. hodgsonii (Flook.f. & Thom.) H.Keng 
in subsect. Blumiana (clade B). The mode of fruit 
dehiscence is a key character to distinguish the two 
subsections; carpels of subsect. Gwillimia show 
longitudinally dehiscence which is widely observed 
in Magnolia, while those of subsect. Blumiana 
show circumscissile dehiscence (Figlar & Nooteboom, 
2004), which means that identification of individuals 
or specimens having no fruits to either subsect. 
Gwillimia or Blumiana is rather difficult, in addition, 
it is hard to obtain good specimens (with fruits) of 
these tree species in tropical forest. Therefore, results 
of our molecular analysis suggest that subsections 
Gwillimia and Blumiana should be abolished and 
only section Gwillimia should be recognized. 
Fruits occurring in Blumiana as well as in South 
American subsect. Talauma probably are a parallel 
adaptation to the tropical rain forest. Apart from 
the fruits the morphology of M. henryi (subsect. 
Gwillimia) and M. hodgsonii (subsect. Blumiana) 
is so similar that, in the herbarium, they are very 
difficult to recognize as separate species (observa¬ 
tion by Nooteboom). More detailed morphological 
studies linked with identification by means of 
molecular techniques may be helpful to solve this 
problem. This may also be true for the other para¬ 
phyletic relationships found in M. liliifera with M. 
villosa (Miq.) K.Keng (clade B). Magnolia liliifera 
is a complex species which contains at least four 


varieties (Frodin & Govaerts, 1996), therefore, our 
finding is supporting the complex status ofM liliifera, 
and we need to re-evaluate the taxonomic statuses 
of the varieties and a related species, M. villosa. 

Clade A was further divided into five clades 
(clades E-I), but the relationships among them was 
not suggested in the MP tree (Fig. 2). The ML tree 
suggested that sister relationship between clades F 
(section Kmeria) and H (section Macrophylla) 
with less than 50% of bootstrap value (Fig. 3). 
Clade E, which was weakly supported in the MP 
and ML trees (60% and 57%, respectively), con¬ 
sists of three well supported subclades (clades J, K, 
and L), which correspond to section Michelia (L.) 
Baill. (clade J; 99% in MP and 100% in ML), and 
sect. Yulania (clade K; 83% in MP and 89% in 
ML), both in subg. Yulania, and subg. Gynopodium 
(clade L; 100% in MP and ML), respectively. 
Section Michelia (clade Ala) taxonomically consists 
of four subsects., Michelia (L.) Figlar & Noot., 
Elmerrillia (Dandy) Figlar & Noot., Maingola 
Figlar & Noot., and Aromadendron Figlar & Noot. 
The molecular phylogenetic trees clearly indicated 
close relationships between subsects. Michelia and 
Elmerrillia (with 91% in MP and 95% in ML), and 
between Maingola and Aromadendron (with 73% 
in MP and 76% in ML), which are also supported 
by morphological features, i.e., the former two 
subsections have flowers on proleptic brachyblasts 
and two to six ovules per carpel, the latter has ter¬ 
minal flowers (rarely proleptic) and two ovules per 
carpel (Figlar, 2000; Figlar & Nooteboom, 2004). 

We included three morphologically putative 
hybrids in the molecular analysis. Because plastid 
(chloroplast) DNA of most flowering plants is 
essentially maternally inherited (Mogensen, 1996), 
the phylogenetic tree presented here may represent 
a maternal lineage of Magnolia. Magnolia x alha 
(DC) Figlar, is a well known taxon and ornamentally 
cultivated in tropical and subtropical Asian coun¬ 
tries. The plant is often sterile and is believed to be 
a hybrid between M. champaca and M. montana 
(Blume) Figlar & Noot. (Frodin & Govaerts, 1996; 
Nooteboom & Chalermglin, 2009), although we do 
not have any molecular or direct evidence to support 
this hypothesis. The phylogenetic tree indicated 
that M. x alha was sister to M. montana (98%). The 
other undescribed putative hybrids, M. champaca 
x haillonii (“Champa khao”; Nooteboom & 


152 


THAI FOREST BULLETIN (BOTANY) 39 


Chalermglin, 2009) and M. champaca x M x alba, 
showed identical sequences with M. champaca, 
indicating mothers of these hybrids should be M. 
champaca, although further evidence is needed to 
further confirm whether they are truly hybrids by 
detailed molecular and morphological studies. 

Phylogenetic relationships among species 
collected from Peninsular Malaysia and Vietnam, 
e.g, M. foveolata (Merr. ex Dandy) Figlar, M. 
mediocris (Dandy) Figlar, M. sirindhorniae, and 
M. scortechinii (King) Figlar & Noot., in subsect. 
Michelia were not suggested. Magnolia foveolata 
and M. mediocris showed the same sequences of all 
regions examined, suggesting they are genetically 
closely related and morphologically recently diver¬ 
sified, but we need an additional examination to 
confirm the sister relationship because we used 
sterile individuals for the samples. Because sect. 
Michelia is the largest group in Magnolia (ca 65 
spp.; ca 50 spp. in subsect. Michelia ; Figlar & 
Nooteboom, 2004), we need more comprehensive 
sampling and sequencing to reveal relationships 
among species in this group and additional taxo¬ 
nomic examination to better understand the species 
boundaries. 

Clade F (100% in MP and ML) corresponds 
to sect. Kmeria (Pierre) Figlar & Noot., including 
M. thailandica from Thailand (Nooteboom & 
Chalermglin, 2009), which is characterized by 
unisexual flowers (only this group). Flowever, a 
sister group of this peculiar group in Magnoliaceae 
was not suggested. 

Clade G (59% in MP and 60% in ML) con¬ 
tains sects. Magnolia (American group; clade M, 
100% in MP and ML), Manglietia (Blume) Baill. 
(100% within clade N) and Rhytidospermum Spach 
(86% and 78% within clade N in MP and NL). In 
the Manglietia clade there are two relatively well- 
supported clades; one clade (86% in MP and 
88% in ML) consists of M. aromatica (Dandy) 
V.S.Kumar, M. grandis (Hu & W.C.Cheng) 
V.S.Kumar, M. hookeri (Cubitt & W.W.Sm.) 
D.C.S.Raju & M.P.Nayar, andM insignis, the other 
(98% in MP and ML) consists of M. conifera 
(Dandy) V.S.Kumar, M. garrettii (Craib) V.S.Kumar, 
M. blumei Prantl. However, relationships among 
the species in each clade are not resolved and mor¬ 
phological characters supporting each clade are 


uncertain. In addition, three individuals (samples) 
of M. utilis (Dandy) V.S.Kumar gave different 
sequences and did not form a single clade, suggest¬ 
ing that we need additional taxonomic and molecular 
works on these species as well as on sect. Gwillimia. 

Clades H (sect. Macrophylla Figlar & Noot.) 
and I (sect. Auriculata Figlar & Noot.) consist of 
temperate American species and had previously 
been treated as section Rhytidospermum sensu 
Dandy (1978), which is characterized by having 
leaves arranged in false whorls. Morphologically 
they are very similar, but molecular data does not 
indicate their close relationship, thus supporting 
their exclusion from subsect. Rhytidospermum sensu 
Figlar & Nooteboom (2004). 

The molecular phylogenetic tree constructed 
in this study provides information on the phyloge¬ 
netic positions of Magnolia species from continental 
Southeast Asia, some of which have not been sug¬ 
gested in previous molecular phylogenetic studies 
(Azuma et ah, 1999; 2001; 2004; Shi et ah, 2000; 
Kim et ah, 2001; Suh & Kim, 2001; Wang et ah, 
2006). In addition, circumscription of taxonomic 
groups (sections and subsections) proposed by 
Figlar & Nooteboom (2004) matches the phylogeny 
(monophyletic clades) except for subsects. 
Gwillimia and Blumiana. However, phylogenetic 
relationships among major clades are still unresolved 
and interspecific relationships within some large 
clades are also not suggested. Additional sequence 
data may improve resolution of the phylogeny, and 
taxonomic re-examination of species boundaries 
noted in this study should be conducted along with 
molecular studies. 
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Table 2. Primer sequences used in this study. 


Foward (5'-3') 


Reverse (5'-3') 


tmK intron 



C7FM 

TGCTAACTCAACGGTAGAGT 

C7RM 

TMF 

TTT CT CTATTAT GGAGT G 

C10RM 

C10FM 

T CCGCTT CTATTT CAGGAGT 

MR2 

MF1 

CT GCT GG ATAC AAG AT GCCC 

MR1 

MF2 

GAAAT GG AA AT CT C ACCTT G 

MR 

MF3 

CGATTT GGGCGGATATAC AG 

TMR 

trnK-39U F» 

GGGGTT GCTA ACT C AACGG 

C7RG1 1 * 

C7FG1 * 2) 

GGT C ACATAGAG A ACGG ATT 

tmK- 2R 2 > 

nclhF 



C11FM 

AGAC ATAT C AATAT GCGT GG 

C11RM 

C11FM2 

TTCCTT GGT GC AGTT GCTAA 

C11RM2 

C11FM3 

C AAAATTAC AGT GGC ACTAAAA 

C11RM3 

C11FG1 3) 

GGGTATAAGAGGATTAGCCG 

350-2R 3 > 

tmH-psbA 



C8F 

CGCATGGTGGATTCACAATC 

C8R 

atpB-rbcL 



C9F 

AGAACCAGAAGTAGTAGGAT 

C9R 

ATM 

CAT GGAT GAATT CT GCCTA 

RBM 

C9F2-138 4 ’ 

GGCTGTT GTTCCTTATTT CAG 

C9R2-138 4 * 



RBM2 5) 


CTAGT CGGAT GGAGTAGATAA 
T CTAGCACACGAAAGT CGAA 
C AAGGT GAG ATTT CCATTT C 
GGGC AT CTT GTAT CCAGCAG 
GTT CGTAAAAAAT CG ATCC A 
C ACT CC ATAATAG AG AAA 
TAGGAAGT GTT GTT GCCG AG 
AACTAGT CGGAT GGAGTAG 

ACATATTTGATACCTTCTCC 
TT CCCC ATATAG ATATT G AA 
T CC AT GGC AT C AGGTAACC A 
GGAAGAAAAGGAGGATCCGG 

AG ACCTAGCT GCTAT CG AAG 

AC ACC AGCTTT G AAT CC AAC 
TAGGC AG A ATT CAT CC AT G 
AATT GACT CTT GAC AGT GAT 
G AAATT GACT CTT GAC AGT GA 


l> for a sequence between trnKy and matK of Galbulimima 

2) for a sequence between matK and trnK3 ’ of Galbulimima 

3) for parts of ndhF sequences of Degeneria and Galbulimima 

4) for M. thailandica instead of ATM or RBM 

5) for M. sirindhorniae and M. scortechinii instead of RBM 
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Figure 1. Relative positions of amplification and sequencing primers on schemed maps of four cpDNA regions. 
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cpDNA MP tree 
(ca. 5800 bp) 
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Figure 2. The strict consensus tree of 334 most parsimonious trees of Magnoliaceae based on trnK intron including matK, trnH— 
psbA, atpB-rbcL, and ndhF sequence data [total of 5960 characters after alignment; tree length = 999; consistency index (Cl) = 
0.8649; Cl excluding uninformative characters = 0.7648; retention index (RI) = 0.8985]. Numbers on branches indicate bootstrap 
percentage (>50%). The taxonomic system proposed by Figlar & Nooteboom (2004) is also presented. 
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cpDNA ML tree 
(ca. 5800 bp) 
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Figure 3. The maxmum likelihood tree of Magnoliaceae based on trnK intron including matK, trnH-psbA, atpB-rbcL, and ndhF 
sequence data (total of 5960 characters after alignment). Numbers on branches indicate bootstrap percentage (>50%). Alphabetic 
marks indicating branches are corresponding to those in Fig. 2. 
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